Exact Results of Dimerization Order Parameter in 
SU(n) Antiferromagnetic Chains 
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Abstract: After a proper definition of the dimerization order parameter for a spin- 
s' system, I show that this order parameter in the SU(n) (n = 2S + 1) antiferro- 
magnetic chains (or equivalently the 577(2) spin-5* chains with Hamiltonians which 
project out singlet states) is, in the thermodynamic limit, directly proportional to the 

3 ■ 

i ■ staggered- magnetization in the corresponding spin-| XX Z chains which had already 

Ch ! been mapped onto the SU(n) chains. 
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The spin-1 antiferromagnetic chain with the pure biquadratic exchange has the 
Hamiltonian 

# = -^T(s,-s l+1 ) 2 , (i) 

i 

where the summation over % runs over all spins with either free ends or the usual 
periodic boundary condition. Parkinson [1] first discussed the possibility of a mapping 
of Eq. (1) onto the spin-i XXZ chain with the anisotropy A = |, which is in general, 
apart from a constant, described by the Hamiltonian 

H = ~ E H ^ H ii+i = ^ + + l^ 1 - Wi). ( 2 ) 

i 

where a a (a = x, y, z) are Pauli matrices and A is the anisotropy parameter. Barber 
and Batchelor [2] later have shown that the Hamiltonian of Eq. (1) with free ends is 
indeed exactly equivalent to the 9-state quantum Potts chain. They then obtained the 
ground-state energy and the excitation gap by the mapping of the Potts chain onto 
the spin-i chains of Eq. (2) with A = |, and with fields ±|\/5 applied to the two 
ends respectively. Kliimper [3] obtained independently these exact results, and he also 
presented results for the correlation length. 

Based on these exact results, Kliimper [3] and Affleck [4] showed that the ground 
state and low-lying excited states of a series of SU(n) (with n = 2S + 1) antiferromag- 
netic chains can all be similarly obtained. In particular, Affleck [4] showed that the 
SU {n) chains with free ends can be mapped, in a similar fashion, onto the correspond- 
ing spin- 1 XXZ chains of Eq. (2). The generic Hamiltonian of these SU(n) chains 
with free ends are given by 

N-l 

# = -E P °+i(^ (3) 

i=i 

where is the number of spins in the chain, and P/j(S) is the projection operator which 
projects out the state with total spin J of the pair and Sj with (Si + Sj) 2 = J( J+l). 
For 5 = Eq. (3) reduces to the usual Heisenberg model plus a constant, whereas for 
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S = 1, Eq. (3) is equivalent to the pure biquadratic chain given by Eq. (1) with free 
ends because Pf i+1 = [(Si ■ S i+1 ) 2 - l]/3. For S = § , Eq. (3) becomes, 

N-l 

H = YY^Y1 t 6 °( s * ' s *+!) 3 + 8 °( s * ' s *+i) 2 " 372S * ' s *+i - 29? ]- ( 4 ) 

i=l 

We note that this form of the Hamiltonian is quite similar to that of the spin-| chain 
proposed by Babujian [5], which is fully integrable by Bethe's ansatz. The mapping 
of Eq. (3) for a general n = 2S + 1 onto the spin-i chains of Eq. (2) is given by the 
relation, n = 2A. In this paper, I shall use the free-end boundary condition as in Ref. 
[2] and [4], and with even total number of spins N in the N — > oo limit. The final 
results are independent of the boundary condition in the thermodynamic limit. 

Although the exact values of the ground-state energy and excitation gap have 
been obtained for the Hamiltonians of Eqs. (1) and (3) for a general n by the mapping 
onto the spin-^ XX Z chains of Eq. (2) which was exactly solved by Bethe's ansatz, and 
although people are convinced [2,4] that the systems described by Eq. (3) are dimerized 
for any n > 2, it does not seem possible to calculate directly the dimerization order 
parameter [6] which is usually defined in the thermodynamic limit, by 

D^USi-i-Si-Si-Si+i)), (5) 

where the angular brackets denote a ground-state expectation. We note that the ab- 
solute value of D is independent of i in the thermodynamic limit (but % should be far 
away from boundaries since the free-end boundary condition is used here). 

In the course of studying spin-lattice dimerization and trimerization problems [7], 
I have come to realize that there is another effective, perhaps more proper, definition 
of the dimerization order parameter for a general spin-S 1 system. This new order 
parameter, which was given by the ground-state expectation of a square matrix with 
dimension given by the number of states for a two-atom spin-S 1 system [7], can in fact 
be equivalently expressed as the ground-state expectation of the projection operator 
if j0 fEq. (3) as, 

D(n) = ((P?_ U (S) - i? m (S))>, n = 2S + 1. (6) 
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For S = this definition of order parameter is identical to the usual definition of 
Eq. (5), but it is certainly different from Eq. (5) for any S > ^, though it is obvious 
that both definitions can effectively measure the order of dimerization. It is clear 
that the key difference between the definitions of Eqs. (5) and (6) lies in the fact that 
operator Si ■ Sj+i of Eq. (5) has in principle a projection of all states with a J value of 
the paired spins with (Si + Si+i) 2 = J(J + 1), unlike Pf i+1 of Eq. (6), which projects 
out only the singlet state of the pair. 

Both Eqs. (5) and (6) for the definition of dimerization order parameter are still 
meaningful for systems with the periodic boundary condition. However, some care 
should be taken since a dimerized system has two degenerate ground states and the 
expectation with respect to an equal admixture of them will yield zero result in Eqs. (5) 
and (6). With the free-end boundary condition, one has the advantage of a non- 
degenerate ground state. 

To see how one comes to the definition of Eq. (6) , we consider the case of perfect 
dimerization for a spin-*? chain. It is convenient to discuss dimerized states in the 
valence-bond basis. Spin operators can be usefully written in terms of two pairs of 
Schwinger bosons as, 

S+ = a+b, S~=ab+, S z = \{a + a-b + b), (7) 

where a, a + and 6, b + obey the usual boson commutation relations. In this representa- 
tion, a spin-S 1 state with S z = m (-S < m < S) is written as 

(a + ) s+m (b + ) s ~ m 
|m)= v^THf^-m)! 1 ^' (8) 

where |0) is the vacuum state of the bosons. A valence bond is simply a spin-singlet 
configuration, which can be written by the so called valence bond operator C^~- , defined 
by 

C+=aib+-ap+- (9) 



For example, the singlet state of two-atom spin-| system is given by a single valence- 
bond configuration, 

l*i> = ^|o> = |i -i>-|-i i> ; (io) 

for S = 1, this singlet state is given by a two-bond configuration, 

l*2> = (C+f\0) = 2(|1, -1) + | - 1, 1) - |0, 0)). (11) 

Generally, the singlet ground state of a two-atom system, each with spin S, is given by 
a 25-bond configuration as 

1*25) = (C+) 2S |0). (12) 
The perfect dimerization state, \D), of the spin- S chain can then be written as 

( i \ f ?4 2 

lD)= { n[(n -I)!]* ) II^-i^lO), n = 2S + l, (13) 

where I have included the normalization factor. It is clear that if the periodic boundary 
condition is used, the dimerization state of Eq. (13) will be doubly degenerate, as 
mentioned earlier. But as we are using the free-end boundary condition here, \D) of 
Eq. (13) is the only choice for the perfect dimer state. The basic algebras in the SU(n) 
chain of Eq. (3) are given by the following two operations, 



^^)(C+) 2S |0) = (C+.) 2S |0); (14) 

and 

p?jw (c+nc^no) = l - (c+.) 2s (c+ ) 2s | ), (is) 

where the four indices k,i,j,l are all different from one another and, as before, n = 
2S + 1. From these two equations, one can in fact prove that the operator, n*Pj ) i+1 (S), 
obeys the Temperley-Lieb algebra [8], which is the key to the mapping of Eq. (3) 
onto the spin-^ XX Z chain of Eq. (2) with the free-end boundary condition [2,4], by 
considering all possible valence-bond configurations of the type of Eq. (12) involving 



four consecutive atoms. Affleck [4] has provided another more elegant proof by using 
the fermion representation. 

One can also define a normalized version of Eq. (6) as 

D'(n) ee i-<(J?_ 14 (S) - P° +1 (S))>, n = 2S + 1, (16) 

where D is the expectation value of (Pi_n(S) — P® i+1 (S)) with respect to the perfect 
dimerized state \D) of Eq. (13). Using Eqs. (13)-(15), it is a straightforward calculation 
to obtain this expectation value as D = 1 — 1/n 2 . 

Since operator n*P- > i+1 (S) obeys the Temperley-Lieb algebra [2,4], one can write 

Pfi+i - llHffli + lsmh9(al +1 - af)], (17) 

where a a (a = x,y,z) are Pauli matrices and Hff?-^ is given by Eq. (2) with A = 
cosh6> = Using this transformation, one can straightforwardly calculate the order 
parameter D'(n) of Eq. (16) by using the mappings of the exact ground-state of Eq. (3) 
onto that of the corresponding spin-| chains of Eq. (2). One thus obtains 

D'(n) = -^V^4(a*), (18) 

where the expectation is with respect to the ground-state of the spin-| XX Z chain 
of Eq. (2) with the anisotropy A = n/2, and where I have used the fact that, in the 
thermodynamic limit, one has 

in the same expectation. Eq. (19) simply reflects the well-known fact that the infinite 
spin- 1 XX Z chain has no dimerization long-range order but a staggered magnetization 
order. Fortunately, this staggered magnetization had already been exactly calculated 
by Baxter [9] twenty years ago as a function of the anisotropy, 

oo 

> /^=l + 2£(-l)»«p(-2n 2 60 

Z 1 (20) 
X>xp[-(n-l/2)V/20], 

n=l 
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where two expressions are equivalent, the first being rapidly convergent at large 9 
while the second at small 9. For S = 1 and §, (of) 0.5028 and 0.7335 respectively 
to the accuracy of four significant figures. Therefore, one has D'(3) ~ 0.4216 and 
D'(4) 0.6776 to the same accuracy of four significant figures for the corresponding 
SU(3) and SU(4) chains of Eq. (3) respectively. One sees also that in the limit of 
S — > oo, (of) = 1 and hence D'(oo) = 1 as expected. 

Based on the definition of the order parameter in Eq. (6), one can define a cor- 
responding 'four-spin' correlation function as ((P^ i+l Pjj +1 )), in similar fashion to the 
usual definition of the four-spin correlation function, ((Si • Sj+i)(Sj • Sj+i)), for the 
order parameter of Eq. (5). Likewise, by taking the long-range limit (i.e., \j — i\ — > oo) 
in ((Pi i+1 Pjj +1 )), one should be able to obtain the value of order parameter D(n) 
(or D'(n)). This is useful if one is to carry out finite-size calculations with periodic 
boundary conditions. 
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